C20: Fullerene, Bowl or Ring? New Results from Coupled-Cluster Calculations 
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Contrary to recent experimental evidence suggesting that the monocyclic ring is the most stable 
20-atom carbon species, highly accurate calculations convincingly predict that the smallest fullerene, 
the dodecahedron C20, has the lowest energy. A related corannulene-like bowl is nearly degenerate 
in energy to the fullerene. Thermodynamic considerations suggest that at formation temperatures 
of around 700 K the bowl should be the dominant species. The recent application of gradient 
corrections to LDA which supported the ring structure is qualitatively in error. 



PACS number: 

The 20-atom carbon cluster is the smallest size that 
can form a closed fullerene molecule. Since small 
fullerenes have been proposed as possible intermediates 
to Ceo or larger fullerenes and because of the rele- 
vance of this molecule to the general question of the sta- 
bility of fullerenes, this cluster has received a great deal 
of attention Jp - [J . Formation experiments by Bowers and 
co-workers g7] and by Jarrold and co-workers || have 
probed the structure of carbon clusters over a wide range 
of sizes, including C2o- In these experiments monocyclic 
rings are dominant for sizes between 11 and 30 atoms. 
There is no evidence for species with a fullerene struc- 
ture for sizes less than 32 atoms. This suggests that 
rings rather than fullerenes are the most stable isomers 
below n=32. However, in these experiments the systems 
are far from equilibrium. It may be difficult to anneal 
the systems to produce the most stable species, therefore 
missing the fullerene. 

Many theoretical calculations of varying sophistication 
have been performed to try to predict stability J^-@,^] . It 
is well known from calculations on smaller clusters, C2- 
C10 that in order to predict the relative energies of carbon 
isomers, very high accuracy methods (few kcal/mol) are 
required. On the other hand, the highest level calcula- 
tions that have been reported on the C20 molecule are 
based on many-body perturbation theory (MP2). Prior 
calculations using these methods for smaller systems have 
shown that for carbon clusters the convergence of pertur- 
bation theory is slow, with reliability of the order of 50 
kcal/mol Isomerization energies in C20 and other 
carbon clusters can be much less than this (roughly 10 
kcal/mol), so higher accuracy methods are necessary to 
make definitive assignments of stability. Methods with 
such accuracy are commonly applied to small molecules, 
but rapidly become computationally intractable as the 
system sizes increases. As we discuss in this letter, they 
are just feasible for the C20 molecule. 



Because of the expense of even low-level calculations 
on systems of this size, attention in previous work and 
in our calculations has focused on three structures which 
appear to be candidates for the most stable configura- 
tion, (Fig. [j]): a monocyclic ring structure (the ring I 
structure dominant in the experiment by von Helden et 
al. 0), a polycyclic corannulene-like bowl comprising five 
fused hexagons which forms part of the Ceo cage, and a 
small fullerene of near dodecahedral structure. Some ear- 
lier estimates of the energies of these structures are sum- 
marized in Table [j]. Lowest level first-principles calcu- 
lations (Hartree-Fock self-consistent field, SCF) Q sup- 
port experimental observations |^| by assigning the ring 
as the lowest energy form. However, perturbational im- 
provements (MP2) Q reverse these predictions, giving 
the fullerene structure the lowest energy. On the other 
hand, local density approximation (LDA) calculations for 
C20 support the fullerene as the lowest energy structure 
p,[§. However, when gradient corrections [|| are made 
to the LDA calculations, the order is reversed, predict- 
ing the ring as the lowest energy geometry in agreement 
with experiments. The change in the ring/fullcrcnc sep- 
aration on going from SCF to the inclusion of correla- 
tion at the MP2 level is roughly 150 kcal/mol. A similar 
large energy difference is present between the gradient 
corrected LDA and uncorrected LDA. Given the uncer- 
tainties inherent in all the computational methods, none 
of the calculations can be considered definitive |fTo|] . We 
note that the reliability of the gradient corrections is a 
question of considerable relevance since this method has 
been suggested as a relatively tractable way to obtain 
results within chemical accuracy (1 kcal/mol). 

We present here results from theoretical calculations 
on C20 at a much higher accuracy than previous work (of 
the order of a few kcal/mol). Since electron correlation 
effects can have a very large influence on the energetics 
of the isomers, we have performed coupled-cluster calcu- 
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lations using the CCSD(T) method JT|]. This method 
includes correlation effects of single and double replace- 
ments of Hartree-Fock orbitals to infinite order, and in 
addition a perturbative estimate of triple replacements. 
Extensive series of calibration calculations have estab- 
lished that the CCSD(T) method is the most accurate 
SCF-based treatment of electron correlation that can fea- 
sibly be applied to a molecule of this size. 

The coupled-cluster method provides an almost exact 
description of the electronic structure of molecules that 
we are interested in here (see below and in Parasuk and 
Almlof jlO|]) provided a sufficiently flexible orbital basis 
is provided. Unfortunately, the increase of the computa- 
tional load with number of orbital functions is so large 
that completeness of the orbital basis is a source of error. 
The calculations here were performed using Dunning 's 
cc-pVDZ basis jl2| on each carbon atom. This is a con- 
tracted set containing three s functions, two p functions, 
and a d function, obtained from a (9s 4p Id) primitive 
set. It provides a "double zeta" description of the va- 
lence electrons (two functions per shell) and a minimal 
description of the Is core electrons, plus a higher an- 
gular momentum function to describe polarization and 
correlation effects. This level of basis set is the smallest 
capable of yielding reliable results in correlated molecular 
calculations. Studies (see e.g., Ref. IllJ) have shown that 
this basis set, optimized for treating correlation, performs 
better than other sets of the same size. Correlation ef- 
fects involving the carbon Is electrons were not included 
because this is expected to have a negligible effect on the 
predicted energetics. 19.5 million parameters were opti- 
mized iteratively in the fullerene CCSD(T) calculation 
and another 29 billion were estimated perturbatively. 
The computational effort is near the limit of the capa- 
bilities of present day supercomputers even for a fixed 
geometry. Therefore, it is not feasible to search for the 
optimum structure at this level of calculation. To define 
input geometries for the CCSD(T) calculations, we have 
determined C20 structures at two levels: SCF and LDA 
|l4|| . The calculations were carried out on the CRAY 
C90 computer at SDSC using the programs molecule- 
SWEDEN Jl5| and titan . Some of the C 6 results 
discussed in this work were obtained with ACES 11 (l7) . 

We should note some points about the symmetry of 
the various C20 isomers. The putative dodecahedral 
fullerene undergoes a Jahn- Teller distortion: previous 
SCF optimization produced a structure with C2 sym- 
metry, whereas the LDA optimization yields a structure 
with Ci symmetry Both optimized structures are close to 
a dodecahedron. The monocyclic ring is found to have 
Cioh symmetry, with alternating bonds and angles, in 
both SCF and LDA optimizations. For computational 
reasons, this was treated within the subgroup Cih in the 
coupled-cluster calculations. The bowl is found to have 
C$ v symmetry at the SCF level, and is very close to this 
symmetry at the LDA level. Only the C s subgroup of 



C$ v was used in the calculations. All calculations were 
done for the singlet spin state: SCF calculations demon- 
strate that all three structures have a closed-shell ground 
state @. 

Energies relative to the fullerene are given in Table |. 
In agreement with previous results 0, at the SCF level, 
the ring is predicted to be the most stable isomer, by 
at least 50 kcal/mol. The bowl is intermediate between 
two other isomers but closer to the ring in energy. The 
situation is quite different using LDA energies, where the 
fullerene is predicted to be most stable, with the bowl 
closer to the fullerene. The MP 2 results are similar to 
those of LDA. Note that the bowl in all these calculations 
is well separated from the fullerene. 

The most important entry in Table | are the CCSD(T) 
results. These results agree with the LDA results, also 
predicting the fullerene to be more stable than the ring. 
However, at this level, the fullerene and bowl are pre- 
dicted to be essentially degenerate at the LDA geometry 
and the bowl is predicted to be lower in energy at the SCF 
geometries. This ordering is sensitive to the remaining 
error in the correlation energy, but the near degeneracy 
of the bowl and the fullerene will remain. 

The CCSD(T) results clearly predict that the fullerene 
and bowl are more stable than the ring. But we still need 
to assess the accuracy of the method. Other calcula- 
tions have shown that for systems for which the Hartree- 
Fock determinent dominates the configuration expansion, 
CCSD(T) recovers essentially all the correlation energy 
available with a given basis set. The scaled norm of the 
single replacements (7i) as calculated by CCSD(T) may 
be used as a guide to the accuracy of a single determi- 
nent description Q. If 71 is less than 0.02 the CCSD(T) 
approach yields results whose reliability depends almost 
entirely on the accuracy of the atomic basis. In these cal- 
culations, T\ was found to be 0.0153, 0.0136, and 0.0151 
for the fullerene, bowl and the ring respectively. Thus, in 
a complete basis, the accuracy will be of the order of a few 
kcal/mol. However, the cc-pVDZ basis used is still quite 
far from complete. For example, it recovers only about 
75% of the valence-shell correlation energy of the carbon 
atom. Presumably, since we are looking for energy dif- 
ferences, much of this error would be shared by the three 
structures and would cancel in the energy difference cal- 
culation, but it is desirable to provide some calibration 
for this. The present calculations have the largest basis 
that we could conceivably use in a CCSD(T) calcula- 
tion for C20, so we can hardly perform the calibration on 
this system. Instead, in order to ascertain the accuracy 
of the results we have performed similar calculations on 
the smaller system C%, which we can treat with higher 
accuracy basis sets. The principal problem is the slow 
convergence of the dynamical correlation energy (that 
treated accurately by CCSD(T)) with angular quantum 
number I of the orbital basis: the behavior is as l~~ A . Our 
cc-pDVZ basis includes up to d atomic orbitals only. In- 
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eluding higher terms such as / and g functions increases 
the computational load very quickly. Nevertheless, for 
the C6 molecule it is possible to explore basis sets up to 
g functions (using again Dunning's correlation-consistent 
sets) |l2| . The results are summarized in Tab le p[ More 
detailed results will be presented elsewhere. ]20| For Cq 
we have examined three structures: a linear triplet state, 
a D^h ring and a D^h ring (lowest energy). The rela- 
tive energies of the two ring isomers do not change with 
increasing basis. However, the relative energy of the lin- 
ear structure changes by a few kcal/mol as the basis is 
increased. 

The first thing to note is that the corrections in Table 
H are of the order of a few kcal/mol. Similar behavior 
with basis set has been observed in C\q. 0] From this 
we can confidently assert that increasing the size of the 
basis used for C20 certainly would be unlikely to change 
the order of the fullerene vs. the ring structures, since 
this is already at least 40 kcal/mol. 

The relative stability of the bowl and the fullerene is 
a more difficult question. At the present level CCSD(T) 
calculation, there is a considerable difference between the 
correlation energy for the different structures (—2.773 Eh 
for the fullerene, —2.699 for the bowl and —2.629 for the 
ring, at the SCF geometries). Since the main error in the 
calculation is in the correlation energy, this implies that 
there is probably a larger energy error in the fullerene 
calculation than in the bowl or the ring. 

Recently, it has been suggested by Martin that the 
angular basis dependence of the energy can be extrapo- 
lated to convergence by a three-term extrapolation for- 
mula based on the broad characteristics of the occupied 
orbitals for a particular molecule. Assuming that the 
majority of the correlation energy available with a given 
basis for a particular system is recovered by a CCSD(T) 
calculation, the additional basis set corrections can be 
estimated by first identifying the types of orbitals, i.e. a, 
tv or lone pairs, and then making corrections which have 
been established for each level of basis set by compari- 
son with accurately measured properties for a test set of 
molecules. This scheme is qualitatively consistent with 
the present results. The results in Table || show that the 
isomerization energy from the D^h to D^h ring changes 
very little with basis, because the orbital character in 
the two molecules is basically the same. On the other 
hand, there are major changes in the number of a and 7r 
bonds in going from the linear to the cyclic forms, giv- 
ing a significant correction in the Martin approach. The 
direction of this calculated correction is consistent with 
the changes observed by explicitly expanding the basis 
sets used (see Table [n]). Unfortunately we cannot cal- 
culate the isomerization energy between the linear and 
ring structures with the very largest basis set, because 
the calculation on the (triplet spin state) linear form is 
not feasible with this large a basis. 

As an additional check on the Martin approach, we can 



refer to the total atomization energy of the C3 molecule 
using the cc-pDVZ basis (see also Ref. Q ) . In this case 
a 30kcal/mol basis set error in the atomization energy is 
predicted by the Martin formula: applying this correction 
brings the predicted atomization energy within the error 
bars of the measured value. 

The stabilities of the C20 isomers estimated via the 
Martin approach are included as the last line in Table Q. 
For both the SCF and the LDA geometries, the fullerene 
is now more stable than the bowl. We note, however, that 
at the SCF geometries the separation is only 5 kcal/mol, 
which is roughly the accuracy of the calculation. There- 
fore, these calculations only weakly support the fullerene 
as the most stable form. 

Another source of error is the lack of CCSD(T) geom- 
etry optimization for each isomer. In the calculations, 
the ring and the bowl had lower CCSD(T) energies at 
the SCF optimized geometries, while the fullerene was 
lower at the LDA geometry. This is naively consistent 
with the greater importance of electron correlation in 
the fullerene, compared to the other isomers. It is diffi- 
cult to say with any certainty what would be the results 
of a CCSD(T) geometry optimization, but if we accept 
the proposition that correlation affects the fullerene more 
than the other isomers, we would expect geometry opti- 
mization to further stabilize the fullerene relative to the 
other isomers. 

The relative stabilities we have presented are all for 
hypothetical vibrationless states at absolute zero. The 
zero-point vibrational energies for the three isomers are 
69.8, 73.3, and 76.8 kcal/mol for ring, bowl, and fullerene, 
respectively, using published SCF frequencies ||. Still 
the fullerene is most stable. However, experimental data 
is taken at high temperature. Therefore, it is desirable 
to explore finite-temperature effects. While there is no 
doubt that these effects will favour the ring isomer due to 
its higher vibrational entropy , the relative energy of the 
ring versus the fullerene is so large that these effects alone 
cannot explain the nonappearance of the fullerene in a 
fully equilibrated system at temperatures less than 1500 
K (estimated from the SCF frequencies). On the other 
hand, even with extrapolations the fullerene and the bowl 
are quite close in energy. Using an average separation of 
9kcal/mol between the bowl and the fullerene, we esti- 
mate that above approximately 700 K the bowl is more 
stable. Below these temperatures the fullerene will be 
progressively more stable. However, it is easy to believe 
that the transition state between any of these structures 
is very high and dynamically difficult to achieve, making 
it very unlikely for the transition from fullerene to bowl 
to occur. This is supported by our dynamical calcula- 
tions which show these molecules to be stable even at 
very high temperatures. 

Finally, our CCSD(T) calculations clearly show that a 
gradient correction to LDA, at least as implemented in 
Raghavachari et al, is incorrect (by 165 kcal/mol!) for 
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this problem. This is unfortunate, since the gradient cor- 
rection approach has been widely advertised as a solution 
to the well-known energetics problems of LDA, and is rel- 
atively easy to implement. We note, however, that the 
uncorrected LDA gives predictions that are semiquanti- 
tatively correct and much better than SCF for a similar 
computational cost. 
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TABLE II. C6 Convergence of isomerization energies with 
basis set size, (kcal/mol) 
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FIG. 1. C20 structures: ring, bowl, and fullerene) 
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